Evolution of the Density of States Gap in a Disordered Superconductor 
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It has only recently been possible to study the superconducting state in the attractive Hubbard Hamiltonian 
via a direct observation of the formation of a gap in the density of states N(u>). Here we determine the 
effect of random chemical potentials on N(u>) and show that at weak coupling, disorder closes the gap 
concurrently with the destruction of superconductivity. At larger, but still intermediate coupling, a pseudo- 
gap in N(ll>) remains even well beyond the point at which off-diagonal long range order vanishes. This 
change in the elementary excitations of the insulating phase corresponds to a crossover between Fermi- 
and Bose-Insulators. These calculations represent the first computation of the density of states in a finite 
dimensional disordered fermion model via the Quantum Monte Carlo and maximum entropy methods. 
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Analysis of the effect of randomness on the super- 
conducting state has a long history, H beginning with 
Anderson's observation Q that in the presence of non- 
magnetic disorder the Cooper wavefunction can be con- 
structed by replacing {k ],—k j) pairs by an appro- 
priate combination of time reversed, but still extended, 
eigenstates. Recent theoretical work || has attempted 
to understand experiments on the superconducting- 
insulator (SC-I) transition in thin films. Here, the super- 
conducting phase is destroyed by strong disorder which 
localizes the electrons entirely. Numerical work || has 
focussed on boson models in which coherence of the pair 
phase is the central issue, and the role of fluctuations 
in the pair amplitude is suppressed. Granular systems 
might be well described by such pre-formed Cooper pairs, 
and, furthermore, scaling arguments suggest universal 
conductivity might find a natural explanation within bo- 
son systems. || 

Despite the successes of numeric studies of interacting, 
disordered bose systems, it is clear that explicit calcu- 
lations for itinerant electron models, in which pairs can 
break apart, are essential. Experimentally, there is ev- 
idence for both SC-I transitions to a Bosonic state ||] 
and to a Fermionic state 0, distinguished by the exis- 
tence or lack of an energy gap, respectively, in the insu- 
lating phase. There is currently considerable controversy 
regarding the meaning of the two insulting phases seen 
experimentally . || Theories implicating both Bosons ||] 
and Fermions ]10[ exist. 

In this paper we describe the evolution of the single- 
particle density of states for the disordered, attractive 
Hubbard Hamiltonian, a model which can interpolate 
between Bose and Fermi limits. We show that in- 
creasing disorder destroys long range pairing correlations 
and drives a superconducting-insulator phase transition. 
However, the density of states N(ui) shows a gap which 
closes with the destruction of superconducting long-range 
order for relatively weak couplings yet retains a gap be- 



yond the critical disorder strength for larger couplings. 
Sweeps across the disorder-interaction phase diagram 
allow us to locate quantitatively the critical coupling 
strengths for gap and pair order formation in a set of 
representative cases. 

The attractive Hubbard Hamiltonian, in the presence 
of random site energies is, 
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Here c\ a is a fermion destruction operator at site i with 
spin a, rij^ = c\ a c\ a , and the chemical potential \x fixes 
the average density (n). The site energies v\ are inde- 
pendent random variables with a uniform distribution 
over [— V/2,V/2}. The interaction has been written in 
particle-hole symmetric form so that fi — corresponds 
to (n) = 1 at all U and T when V = 0. The lattice 
sum (ij) is over near neighbor sites on a two dimensional 
square lattice. 

We solve for the equilibrium properties of this Hamilto- 
nian using the determinant QMC method. |ll]] Previous 
numerical studies of the clean model have determined the 
phase diagram jl2| by a finite size scaling analysis of the 
equal time pair-pair correlation functions, 



P s (j) = i£(A s+J At>, 
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The phase diagram consists of a state with simultane- 
ous charge density wave (CDW) and superconducting 
correlations at T = and half-filling, and a finite tem- 
perature Kosterlitz-Thouless transition (with a maximal 
T c k O.lt for \U\ = 4t) to a phase with power law de- 
cay of pairing correlations off half-filling. In this paper 
we will study fillings (n) — 0.875 for which strong CDW 
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correlations are absent and the transition temperature to 
the SC phase is nearly maximal. 

In Fig. 1 we show the suppression of long range correla- 
tions in P s (j) with increasing disorder strength. A finite 
size scaling analysis of the structure factor determines the 
critical value V c — (3.25 ± 0.2)t for the destruction of the 
superconducting state. V c can also be identified by the 
superfluid density D s and the conductivity. These 
quantities give consistent values, V c = (3.2 ± 0.7)t and 
V c = (3.5 ± 0.5)i, respectively, for the critical disorder. 

lO 
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FIG. 1. The equal time pair correlation P s (j) is shown at 
T = 0.10 as a function of separation j for different degrees of 
disorder. P s {j) remains finite at large j for V = 2.0, but goes 
to zero rapidly with j at V = 4.0. 

The density of states is a quantity which, like the con- 
ductivity, is directly accessible experimentally. Its evalu- 
ation requires an inversion of the integral relation, 

g (-) = ^B c pW4(o)) = /7^t+-^' (3) 

between the density of states N(u>) and the imaginary 
time dependent Greens function G(r). We do this using 
Bryan's method [[L4| and Classic Maximum Entropy [ 15[ , 
using the full imaginary-time covariance matrix. 

Although this approach is fairly well-developed for 
clean systems, the inclusion of disorder raises a number 
of new questions of principle. The proper treatment of 
errors and correlations in the QMC data was the central 
achievement of the ME approach, [|l6| and it is not obvi- 
ous how additional fluctuations from disorder averaging 
will affect the procedure. Therefore, we have checked 
our calculation against analytic results both at weak and 
strong coupling. Besides verifying the numerics, the be- 



havior in these two limits also presages the two distinct 
types of insulating behavior we see in the full model. 

In the non-interacting limit we can diagonalize the 
Hamiltonian and obtain E a , the single particle eigenen- 
ergies and thence N{tu) — J2 a ™ E a ). In this limit 
the determinant QMC method computes the imaginary 
time fermion Greens function exactly. Tests of the ME 
approach in the clean case require that some noise be 
added to the exact G(r) coming from the QMC. g7| In 
the random case, this ad hoc noise does not need to be 
added, since G(t) already has error bars from disorder 
averaging. We find that the ME approach exactly tracks 
the disorder induced broadening of N(ui). Jll| For small 
\U\ we expect disorder to broaden and eliminate the su- 
perconducting gap as well. 

The strong coupling, t = 0, limit also can be solved 
analytically. For a single site, 

e -/3(7/4 e -rA_ + e f(U/4+v i ) e -rA + 

G(T) = 2(e/ 3c/ / 4 cosh/3t> i + e-W 4 ) ' 

(1 + e-P A -)e-P u / 4 
N{LJ) = 2(e^/ 4 cosh^ i + e-^/ 4 ) (5( " " A " } 

H + e -f3A + \ e f3(U/4+Vi) 

+ 2( e W 4 cosh/^ + e-W 4 / (w ~ A+) - (4) 

Here A± = ±U/2 + Vi. After disorder averaging, the two 
delta functions in the density of states N(ui) are broad- 
ened to two distributions centered about ±£7/2, each with 
width V. When V — U these merge, and the gap in N(ui) 
at u> — is closed. However, the thermal factors greatly 
suppress N(u>) near u> = 0. At t = there is no long range 
pairing order, so that this "pseudo-gap" in the density of 
states reflects the tendency for on-site singlet formation. 
We have shown that our QMC+ME procedure accurately 
reproduces the analytic result Eq. 4 as a function of tem- 
perature, disorder strength, and interaction strength. JlJ] 
A central conclusion of our work is that this pseudogap 
behavior persists far from the strong-coupling limit as 
the hopping t is turned on. 

We now study the cross-over between these two possi- 
ble types of effect of disorder on N(oj) in the full model. 
In Fig. 2 we show the evolution of N(u>) with increasing 
interaction strength \U\ at fixed disorder V = 2t. N(ui) 
evolves from its gapless, noninteracting form to possess- 
ing a well-formed gap at \U\ = At. Similarly, in Fig. 3 we 
show the evolution of N(u>) with increasing disorder V, 
at fixed values of the interaction, \U\ — 4t. In contrast 
to the non-interacting case where increasing V results in 
the expected broadening of N(uj), here the gap in N(uj) is 
remarkably robust. Note the values of V are well beyond 
the point where the superconducting-insulator transition 
has occurred, as indicated by the results in Fig. 1 and re- 
lated measurements. [ [l3"| In the repulsive Hubbard model 
at the same value of U/t the gap in N(u>) requires long 
range antiferromagnetic order be present, specifically, the 
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correlation length £ a f must exceed the lattice size. Q In 
contrast, here we find that a gap in N(w) persists despite 
the clear destruction of such long range order. 
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FIG. 2. The density of states N(u) at t = 1, T = 0.125, 
V = 2, |t/| = 1,2, 3, 4 on an 8 x 8 lattice. For sufficiently 
small \U\, scattering by the disorder inhibits phase coherence 
and the superconducting gap is destroyed. A superconducting 
gap forms with increasing interaction strength U . A study of 
the pair correlations P s (j) indicates that gap formation occurs 
precisely with the onset of superconductivity. 

How can we interpret these results? In the weak- 
coupling limit, disorder broadens N(uS) at all energies. 
In particular, the superconducting gap in N(oj) is closed 
at the same point where disorder destroys the long-range 
pair correlations. This closing of the gap is indicative of 
the predominance of fermionic elementary excitations in 
the insulating phase. In the strong coupling limit, the 
density of states is not broadened by disorder near ui = 
and a pseudogap survives. The results of Fig. 3 imply 
that already at a value of interaction strength \U\ = At 
which is only one-half the bandwidth, W = 8t, the super- 
conducting gap is robust to disorder. In previous work, 
p0| it was found that at V = and \U\ — At the spin 
susceptibility exhibits a reduction at low temperatures, 
indicative of bose-like behavior, despite the fact that the 
temperature dependence of the chemical potential fi(T) 
is still clearly that of a degenerate Fermi system. Wc 
conclude here that for an interaction strength as low as 
\U\ = At, the density of states reflects bose-like charac- 
ter at all V. As shown in Fig. 3, at either side of the 
critical disorder V c for the destruction of long-range or- 
der in this coupling regime, there are no indications of 
fermionic excitations in N(u>). 
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FIG. 3. The density of states N(lo) at t = 1, U = -4, 
T = 0.125, V = 0.0, 2.0, 4.0 on an 8 x 8 lattice. A study of 
the equal time pair correlation shows that for V < V c ~ 3.5 
the system is a superconductor. For V > V c the system is 
an insulator. The gap in N(ui) persists even after long-range 
pair correlations, and superconductivity, are destroyed. Fur- 
thermore, there is no enhancement in N(ui) as one approaches 
V c from either side, indicating that the system is not a Fermi 
liquid at the SC-BI transition. 

The phase diagram is shown in Fig. 4. For weak inter- 
actions, there is a SC-Fermi Insulator (SC-FI) transition 
with disorder as seen in one class of experiments. For 
a stronger, but nevertheless intermediate coupling, the 
phase diagram already exhibits the SC-Bose Insulator 
(SC-BI) transition seen in another class of experiments. 
H From Eq. 4 it follows that for small t the pseudo- 
gap in N(w) opens even for very small \U\/V at large 
(3, indicating that the FI-BI crossover has the trajectory 
shown. []l8f Theoretically, a universal resistivity is pre- 
dicted by a bosonic treatment of the transition. |J It 
has been suggested that the presence of fermions causes 
non-universal properties at the SC-I transition. It is 
possible that the attractive Hubbard Model will exhibit 
universal resistivity only in the Bose-Insulator regime, 
though considerable further work needs to be done to 
settle this issue. 

The interplay of finite-size effects and the development 
of a gap in N(ui) is a subtle issue. As described above, 
in the case of the half-filled repulsive Hubbard model 
at U = At, an antiferromagnetic gap in N(ui) at finite 
temperature disappears when the lattice size is increased 
beyond the correlation length of antiferromagnetic order. 
We have studied this issue carefully by further evaluation 
of N(uj) on a range of lattices, and conclude that the gap 
we observe is a robust feature and not an artifact of finite 
lattice size. A crucial difference here from the case of the 
half-filled repulsive model is that the antiferromagnetic 
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transition takes place only at T = 0, whereas the su- 
perconducting transition off half-filling in the attractive 
model takes place at finite temperature. 
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FIG. 4. Phase diagram as a function of disorder V/t and 
interaction strength \U\/t. The solid circles indicate the lo- 
cation of the entry into the superconducting phase as deter- 
mined by our simulations. The remainder of the solid line 
indicates an estimate of the location of the boundary of the 
superconducting region. The arrows correspond to Figs. 2 
and 3, open triangles indicate gapped regions while closed 
triangles indicate no gap in N(u). In the superconducting 
regime the pair correlations have long range order and there 
is a gap in N(uj). In the insulating phases there is no pair 
long range order. The hatched region demarks schematically 
the cross-over between the Bose and Fermi Insulators. 

In conclusion, we have computed the density of states 
N(lu) in a disordered interacting fermion model with 
a combination of the QMC and ME methods. The 
disorder-induced broadening expected in the weak- 
coupling BCS limit is already absent by the time \ U\ = 4t, 
half the single particle bandwidth W = 8t. Disorder 
closes the gap in the density of states for weak interac- 
tions but for stronger interactions has little effect on the 
gap well beyond the point where pair coherence vanishes 
and the system is insulating. At the coupling where a 
metallic state has been previously identified, the density 
of states suggests bosonic excitations on either side of the 
SC-I transition. 
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